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ABSTRACT 


Structural  acoustics  often  utilizes  Fourier 
transformation  to  either  reveal  a  phenomenon  that  is  more 
directly  recognized  in  one  domain  than  in  its  Fourier 
conjugate  domain  or  to  employ  complementarity  in  order  to 
decipher  a  phenomenon.  In  this  report,  a  novel  computational 
technique  is  introduced.  This  technique  advantageously 
performs  this  transformation  in  situations  in  which  specific 
spatial  scales  govern  aliased  factors  in  the  integrand  of  a 
Fourier  transformation.  The  technique  can  factorially  save 
computational  steps  over  corresponding  computations  that 
employ  discrete  Fourier  transform  procedures. 


I.  INTRODUCTION 

Structural  acoustics  deals  largely  with  the  description  of  wave  propagations  in 
structures  and  in  the  fluids  that  are  in  contact  with  their  internal  and  external  surfaces. 
Models  of  the  structure  are  constructed  so  that  they  facilitate  the  derivation  of  these 
descriptions.  A  description  may  be  more  readily  derived  when  the  model  is  idealized  and 
adapted  to  accommodate  a  given  Fourier  domain.  Yet  a  phenomenon  in  the  description 
may  be,  in  some  cases,  more  advantageously  interpreted  in  a  domain  that  is  a  Fourier 
conjugate  to  that  in  which  the  description  is  originally  derived.  Moreover,  a  comple¬ 
mentary  description  of  this  kind  is  itself  of  considerable  benefit  in  the  interpretation  of  a 
phenomenon  in  the  description.  In  this  sense,  one  is  often  called  upon  to  perform  Fourier 
transformations  in  structural  acoustics.  The  purpose  of  this  report  is  to  propose  a  novel 
computational  technique  to  accomplish  a  certain  class  of  Fourier  transformations.  It 
transpires  that  a  number  of  phenomena  in  structural  acoustics  are  featured  in  this  class;  in 
particular,  those  phenomena  that  are  related  to  aliasing  and  to  pass  and  stop  bands  in 
regularly  ribbed  structures  [1-11]. 
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A  typical  Fourier  transform  in  this  report  is  expressed 


fix,  CO) 

• 

j*  dk  exp{-ixk)  F(k,co)  , 

(la) 

F(k,(0) 

• 

f  dx  expiixk)  f(x,co)  , 

(lb) 

where  the  -domain  and  the  x  -domain  are  Fourier  conjugate  domains;  (k)  and  (jc)  are 
a  Fourier  variable  pair;  and  {(o)  is  a  silent  variable  in  these  transformations.  For  these 
equations  to  be  valid,  F(k,co)  and  f(x,co)  need  to  be  well  behaved  functions  of  (k)  and 
(jc),  respectively,  so  that  convergence  can  be  guaranteed.  In  this  report  the  concern  is 
mainly  focused  on  Eq.  (la).  Using  the  self-aliased  wavenumber  operator  Sf^ik),  Eq.  (la) 
may  be  cast  in  the  form 


f(x,co)  =(271)-^'^  [  ^  dk  Si,(k)  {exp(-ixk)F(k,co)]  ; 

O 

KTi  =  (27t)m  .  (2a) 

The  self-aliased  wavenumber  operator  5^  (A:)  is  defined 

Shik)  {N(k,(0)}  =  N(k  +  K  CO)  =  M(k,co)  =  M(k  +  Kj,(o)  ; 

Si,ik)  {K(k,(0)  G{k,co)}  =  K{k,(o)  S^ik)  {G(k,co)}  ,  if 

K(k,co)  =  K(k  +  Kj,co)  ,  (3) 
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where  {b)  is  an  appropriately  selected  spatial  scale  factor, 


=  <1^1  ;  =  j^\  ;  =  i27t lb)  ,  (4) 

and  N(k,(o),  K(k,co)y  and  G(k,Q))  are  arbitrary  but  well  behaved  functions  of  the 
wavenumber  (k).  The  quantities  K(k,co),  by  definition,  and  M(k,(o),  by  construction, 
are  aliased;  i.e.,  once  K{k,co)  and  M{k,(o)  are  specified  in  a  segment  that  spans  the 
width  (k^)  in  the  ^-domain,  they  are  specified,  by  repetitive  extrapolation,  throughout 
that  domain.  In  particular,  if  K(k,co)  and  M(k,co)  are  known  within  the  segment 
0  <  k  <  they  are  known  throughout  the  k  -domain.  Equation  (la)  is  transformed 
into  Eq.  (2a)  by  constructing  an  integrand  that  is  “aliased  in  (k)  with  respect  to  (/Ci)” 
and,  therefore,  the  integral  over  a  single  segment  (kTj)  wide  suffices;  all  the  information 
that  is  in  [exp(-ixk)  F(k,(o)}  over  the  entire  ^-domain  is  concentrated  in 
■S'fc(^)  [exp(-ixk)  F(k,co)}  over  the  segment  0  <  ^  <  of  that  domain.  After  the 
aliased  integrand  is  sampled  in  this  segment,  the  integration  over  this  segment  may  be 
evaluated  adaptively  or  otherwise.  Alternatively,  the  segmental  span  (k^  )  may  be  chosen 
to  be  (k[)  which  is  small  enough  so  that  a  “mean  value  theorem”  can  be  relied  upon  to 
evaluate  the  integral  in  Eq.  (2a).  At  this  stage,  the  alternative  procedure  is  selected  for 
the  evaluation  of  Eq.  (2a);  namely 

f(x,co)  =  (k{)  Sy(k^)  {exp(-ixk^)  F(k^,co)]  ; 

k{  =  {271  lb')  ;  k„  =  (aK{)  ; 

{N{k,co)]  =  T^{N  +  k'  CO)  ;  k'  =  qk'^  ,  (2b) 
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where  (a)  is  a  positive  constant  less  than  unity,  [cf.  Eqs.  (3)  and  (4).]  To  define  “small 
enough,”  the  details  in  the  integrand  in  Eq.  (la)  may  be  conveniently  typified  by  a 
wavenumber  (ky )  and  a  loss  factor  (rjy ) .  The  wavenumber  (ky )  locates  the  position  of 
the  (y)th  detail  and  (jjy  )  estimates  the  extent  of  that  detail  in  the  k  -domain.  In  order  to 
capture  all  relevant  details  in  this  integrand,  when  Eq.  (la)  is  converted  into  Eq.  (2b),  it  is 
necessary  to  impose 

K{  «  (7]yky)  ,  (5) 

for  all  relevant  y’s.  Indeed,  under  the  condition  stated  in  Eq.  (5),  and  after  carrying  out 
the  operation  dictated  by  Eq.  (2b),  (k^)  can  be  safely  set  equal  to  zero  by 

choosing  (a)  to  be  zero.  In  this  limit,  Eq.  (2b)  becomes  a  Fourier  series  expansion  for 
f(x,(o).  Clearly,  the  infinite  summation  that  is  dictated  in  Eq.  (2b)  by  the  operator 
Sb’iko)  cannot  be  met  in  practice;  the  summation  needs  to  be  truncated.  To  ensure  that 
the  tmncated  terms  can  be  discarded,  a  maximum  wavenumber  (kj^ )  is  defined  so  that 
the  convergence  of  F{k,(o)  is  substantially  completed  when  this  wavenumber  is  reached. 
A  number  (/')  is  then  defined  and  it  is  required  that 

f(x,(0)  (*•;)  si;.(k,)  [exp(-ak,)  F(k^,(0))  ; 

kn  =  U'K[)  .  (2c) 

constitutes  an  adequate  approximation  to  Eq.  (2b),  where 


sHik)  [N(k,<o)]  i;  m  +  K'j,ay.  k'j  =  .  (6) 

=  -r 
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[cf.  Eqs.  (2b)  and  (3).]  How  can  an  adequate  approximation  be  defined  in  the  context  of 
Eq.  (2c)?  For  this  purpose  one  needs  to  specify  the  expected  “overall  spatial  span”  {Xf^) 
and  the  “spatial  sampling  span”  (jcf)  for  the  Fourier  transform  f{x,(o).  The  specification 
of  these  spatial  parameters  defines  a  number  (/' )  given  by 

~  ^i)  •  (7a) 

The  uncertainty  principle  and  sampling  theory  relate  }  to  {jc^ ,  }  in  the  form 

=(^mO  ^  2;r  ;  J'^  =  J'  >  27t  {x{k{)~^  .  (7b) 

Summarizing  Eqs.  (2)  through  (7)  yields 

f(x,co)  (/cO  Si(,ky)  {exp(-ixk^)  F(k^,(o)}  ; 

ko  =  (OCK{)  ,  (8) 

and  the  approximation  is  adequate,  provided 

(Vy)''  «{kylK{)  «  J'-,  J’  ={kulK[)  =(Xulx[),  (9) 

is  satisfied,  where  {ky )  and  (rjy )  define  all  details  in  F{k,(o)  in  the  k  -domain  that  need 
to  be  captured  to  render  f{x,(o)  a  reliable  Fourier  transform.  Usually  {x[  ?£•()  is  selected 
as  small  compared  with  unity  and  («)  is  selected  as  zero.  Choosing  the  equality  sign  in 
Eq.  (7b)  and,  as  such,  substituting  it  in  Eq.  (8),  the  resulting  equation  can  be  readily 
manipulated  to  conform  with  the  discrete  Fourier  transform.  The  discrete  Fourier 
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transform  is  the  basis  for  the  formulation  of  the  Fast  Fourier  Transform  (FFT) 
computational  procedures  [12].  In  this  report,  however,  a  different  computational  path  is 
selected. 

II.  FORMALISM  OF  A  DOUBLE-SUM  COMPUTATIONAL  PROCEDURE 

The  approximation  that  underlies  Eq.  (8)  requires  Eq.  (9)  to  be  satisfied.  The 
latter  equation  implies  that  some  knowledge  of  F(k,(o)  is  on  hand  in  that  the  various 
vectors  that  typify  details  in  this  function  can  be  reasonably  estimated.  These 

estimations  are  then  used  to  determine  the  vector  [k^,/'}  that  needs  to  be  imposed  on 
Eq.  (8).  This  vector,  in  turn,  determines  the  “wavenumber  sampling  span”  (k[)  and  the 
“overall  wavenumber  span’  =  The  first  determines  the  separation  between 

adjacent  sampling  positions  and  the  second  ensures  convergence.  Situations  may  arise  in 
which  the  wavenumber  sampling  span  is  selected,  for  some  reason,  to  be  (Kj)  instead  of 
(k[)  and  (k’i)  >  (fc{).  Were  this  selection  of  (Kj)  for  (k[)  to  be  made,  what 
modification  needs  to  be  instituted  in  Eq.  (la)  to  bring  in  a  finite  summation,  while 
ensuring  that  an  adequate  approximation  is  installed?  Indeed,  the  approximation  is 
required  to  match  the  accuracy  that  is  enjoyed  by  Eq.  (8)  when  this  equation  is  covered 
by  the  imposition  of  Eq.  (9). 

Using  the  self-aliased  operator  Eq.  (la)  may  be  recast  in  the  form 

/(x,to)  =  (2;r)~^'^  dJt  Si,(k)  [exp(-ixk)  K(k,co)  G(k,(o)}  ; 

^  O 

=  (271  lb)  ,  (10a) 
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where,  for  the  sake  of  flexibility,  F(k,(0)  is  expressed  in  the  factorial  form 

F(k,co)  =  K{k,(o)  G(k,co)  .  (11) 

[cf.  Eq.  (2a).]  Ensuring  the  convergence  of  F(^,co)  by  an  overall  wavenumber  span 
ki^  =  /Kj,  the  infinite  summation  implied  by  Sf^ik)  is  then  replaced  by  a  finite 
summation;  namely,  Eq.  (10a)  is  approximated  in  the  form 

f{x,co)  dk  slik)  [exp{-ixk)  K{k,(o)  G{k,(o)}  ; 

w  O 


(10b) 


where 


j 

sl{k)  (JV(<:,<B))  =  2  N{k  +  K^,(oy,  Kj  =  jKi  ; 

j  =  -J 


=  ilKlb)  ,  (12a) 

[cf.  Eqs.  (2c)  and  (6).]  With  the  help  of  Eq.  (7)  it  emerges  that  the  approximation  in 
Eq.  (10b)  matches  that  in  Eq.  (2c)  provided  the  equalities  and  the  inequality 

=  A  ;  ;  Icm  = 

J  >  27r(XiKi)-^  ,  (13) 
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are  maintained.  Since  the  new  wavenumber  sampling  span  (k^)  is  larger  than  the  old 
wavenumber  sampling  span  (k{),  the  integral  cannot  be  evaluated  by  a  single  segment  in 
the  manner  of  Eq,  (2c).  Each  segment  in  the  summation  needs  to  be  spanned  more  finely. 
In  this  vein,  the  integral  is  carried  out  over  (1  +  R)  regular  intervals,  notwithstanding  that 
an  adaptive  procedure  may,  under  certain  circumstances,  be  preferred.  Executing  the 
integral  in  the  manner  just  specified,  one  derives  from  Eq.  (10b)  the  Fourier  transform  in 
the  form 


R 

f(x,co)  (l  +  K)"'  (K-i)  2  sl(k,)  {exp(-ixk,)  K(k,.(0)  G(K,(o)}  ; 

r  =  o 

;  K  =  (1  +  ^)"^  =  (rK^)  ;  KTi  =  {2n/b)  . 

(14) 

The  “effective”  wavenumber  sampling  span  in  Eq.  (14)  is  then  given  by  +  k^. 

Using  Eqs.  (9)  and  (12),  the  approximation  in  Eq.  (14)  adequately  matches  that  of  Eq.  (8) 
provided 

{-n^il  +  R)}-^  «  (ky/K^)  «  J  ;  J(l  +  R)  cr  J'  ,  (15) 

notwithstanding  that  the  number  of  terms  in  the  double-sum  in  Eq.  (14)  is 
{(2/  +  1)  -I-  (l-h/?)}  versus  (2/' -hi)  in  the  single-sum  in  Eq.  (8).  However,  the  number 
of  computational  steps  is  {(27  +  1)  (1  +  /?)}  in  Eq.  (14)  versus  (27' +  1)  in  Eq.  (8),  which 
is  in  agreement  with  the  second  of  Eq.  (15)  when  7  and  7'  »  1.  There  is  then  no 
factorial  saving  in  computational  steps  between  Eq.  (14),  as  it  stands,  and  Eq.  (8);  using 
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one  or  the  other  equation  is  merely  a  matter  of  convenience.  The  crux  of  the  report  is 
revealed  when  one  or  the  other  factors  composing  K{k,(o),  etc,  is  aliased  in  {k) 

with  respect  to  (k^),  and  (k^)  matches  the  (Kj)  in  Eq.  (10a).  Under  this  condition,  a 
factorial  saving  of  computational  steps  may  be  accrued  by  an  appropriate  manipulation  of 
the  integrand.  A  few  examples  in  this  regard  are  considered  next. 

m.  INTEGRAND  COMPOSED  OF  ALIASED  FACTORS 

The  expression  depicted  in  Eq.  (11)  is  specialized  in  that  the  factor  K{k,co)  is 
assumed  to  be  aliased  in  (^)  with  respect  to  (k^).  [cf.  Eq.  (3).]  Using  the  properties  of 
the  self- aliased  operator  Sf,{k)  as  defined  in  Eq.  (3),  the  integrand  in  Eq.  (10a)  can  be 
equivalently  expressed  in  form 


Sjjik)  [exp{-ixk)  K(fc,£o)  Gik,(o)]  =  K{k,co)  S^(^)  {exp(-ixk)  G(k,co)}  ,  if 

K(k,co)  =  K(k  +  Kj,(o)  .  (16a) 

In  terms  of  the  approximation  stated  in  Eq.  (14),  the  “R-summand”  equivalence  assumes 
the  form 

Slik^)  [exp{-ixk^)  K(k^,co)  Gik^^co)]  =»  K(k^,co)  {exp(-ixk^)  Gik^yCO)} 

~  ’  (16b) 
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where  kf^-^  is  the  overall  wavenumber  span  with  respect  to  the  convergence  of  the 

j 

operand  under  )  in  this  equation  and  in  general 


{N(.k„(0)} 


N(kr  +  Kj ,  CO)  . 


(12b) 


It  is  noted  that  if  the  right  side  of  Eq.  (16b)  is  used  in  Eq,  (14)  instead  of  the  left,  a 
factorial  saving  equal  to  the  smaller  of  + 1)  and  (2/  + 1)  of  computational  steps  is 
accrued.  Usually,  however,  =  J  and,  therefore,  the  two  factors  are  equal. 
Moreover,  if  the  aliasing  in  K(k,  co)  can  be  cast  in  the  form 


K(k,(0)  =  S,(k)  {KA{k,(0)  Ks(k,co)]  , 


(17a) 


then  in  terms  of  the  approximation  stated  in  Eq.(14),  the  R-summand  equivalence 
becomes 


{exp(-ixkr)  K(kr,(o)  G(kr,(o)} 

=*•  (^,1  (*,,<»)  Kg(k,,0)))  sl'(k,)  [expi-ixk,)  G(Jc,,(0))  ; 

km  =  UiKi)  ,  (17b) 

where  kf^2  the  overall  wavenumber  span  with  respect  to  the  convergence  of  the 
operand  {K/^{k,co)  Kg{k,(o)}  under  S^^(k).  It  is  noted  that  the  factorial  saving  of 
computational  steps  pertaining  to  Eq.  (17b)  remains  unchanged  from  that  derived  from 
Eq.  (16b).  If,  in  addition  to  the  form  of  K(ky  co),  as  stated  in  Eq.  (17a),  it  is  assumed  that 
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K^{k,(o)  is  aliased  in  {k)  with  respect  to  (Kj),  then  the  factor  K(k,co)  can  be  factorized 
further  in  the  form 

K(k,co)  =  K^{k,co)  S,(k)  {Ks(k,co)}  ,  if 

K^ik^co)  =  K^(k  +  Kj,o))  .  (18a) 

In  terms  of  the  approximation  stated  in  Eq.  (14),  the  R-summand  equivalence  assumes 
the  additional  factorization,  over  Eq.  (17b),  of  the  form 

si(k,)  {exp(-ixk,)  K{k^,(o)  G(k,,co)} 

=>  K^{k^,co)  sl^ik^)  KQ{_k^,co)}  [exp{-ixk^)  G{k^,(o)}  ; 

^M3  =  (-^s^i)  .  (18b) 

where  is  the  overall  wavenumber  span  with  respect  to  the  convergence  of  the 
quantity  K^(k,(o)  in  this  equation.  It  is  noted  that  if  in  Eq.  (14)  the  right  side  of 
Eq.  (18b)  is  used  instead  of  that  of  Eq.  (16b),  a  factorial  saving  equal  to  the  smaller  of 
(2^2  +1)  and  (2/3  +1)  of  computational  steps  is  accrued.  Usually,  however,  J ^  =  J 
and,  therefore,  the  two  factors  are  equal.  Under  this  equality,  it  is  noted  that  if  the  right 
side  of  Eq.  (18b)  is  used  in  Eq.  (14)  instead  of  the  left,  a  factorial  saving  of  (27  +  1)^  of 
computational  steps  is  accrued.  In  the  event  that  the  aliasing  in  the  factor  K^{k,(o),  as 
stated  in  Eq.  (18a),  stems  from  its  functional  dependence  on  an  aliased  quantity;  e.g. 

Kj,(,k,(0)  =  /5:j5,(/:)  [Kj^iKco)]]  ,  (19a) 
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the  approximation  stated  in  Eq.  (14)  supports  the  R-summand  equivalence  of  the  form 


Slik,)  {exp(-ixk^)  K(k,,(o)  G(k,,o})} 

=>  K}(k„o>)  Sj’a,)  {Ks(,k„0))]  sIhK)  G(k,.(0))  ; 

KiHKM  =  K^isi‘(.k,)  {K^(k„iom-.  *«4  =  (Aif.).  (19b) 

where  kj^f^  is  the  overall  wavenumber  span  with  respect  to  the  convergence  of  the 
quantity  K^{k,co)  in  this  equation.  It  is  noted  that  there  is  no  factorial  saving  of 
computational  steps  between  the  use  of  Eqs.  (18b)  and  (19b).  If  further  factorization 
exists,  either  in  K^(k,co),  in  Kg(k,co)  or  in  both,  in  the  format  just  prescribed,  one  may 
gain  more  and  more  factorial  savings.  For  the  purpose  of  this  report,  however,  the 
factorization  and  dependence  considered  in  Eqs.  (18a)  and  (19a)  suffice. 

A  special  case  of  Eq.  (11)  is  that  in  which,  in  addition  to  the  aliasing  of  K{k,(o)  in 
ik)  with  respect  to  (Kj),  the  function  G(J(.,(0),  in  this  equation,  is  independent  of  the 
wavenumber  {k)\  namely 

F{k,(o)  =  G^{(o)  K(k,co)  ;  K(k,Q))  =  K(k  +  Kj,co)  .  (20) 

Substituting  Eq.  (20)  in  Eq.  (10a)  derives 

fix, CO)  =  St,ix)  {5(x)  f,ix,Q})]  ; 


Mx,co)  =  [G,ico)/K,]i27r)-^'^  f  '  dk  exp(-ixk)  Kik,co)  ,  (21) 
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where  Si,(x)  is  a  self- aliased  spatial  (x)  operator  with  respect  to  the  fundamental 
separation  distance  (ft);  ft  =  {2k  Ik^).  This  operator  is  defined 

Sf,{x)  {n{x,(o)]  =  S  n{x  +  X:,(o)  ;  jc.-  =  (yft)  ;  =  0  .(22) 

j 

[cf.  Eq.  (3).]  Executing  the  integral  in  Eq.  (21)  in  a  manner  analogous  to  that  resulting  in 
Eq.  (14),  one  obtains 

R 

fs{x,co)  =  {iKy^'^  [G^{q}) /(1  +  R)]  {exp{-ixky  K{k,,co)}  ,  (23) 

r  =  c> 

and  the  evaluation  of  y{x,ct))  is  considered,  in  view  of  Eq.  (21),  to  be  synonymous  with 
the  evaluation  of  f{x,co),  notwithstanding  the  acknowledgment  that,  unlike  f{x,co),  the 
evaluation  of  fj.{x,co)  is  meaningful  only  for  values  of  x  =  x^;  Xj  =  (jb).  In  this 
view,  the  recognition  of  the  aliasing  in  K{k,(o),  as  stated  in  Eq.  (20),  is  worth  a  factorial 
saving  of  {2J  + 1)  of  computational  steps.  This  is  similar  to  the  factorial  saving  that  is 
accrued  in  Eq.  (16b),  utilizing  the  same  recognition.  Indeed,  the  factorization  of  K{k,  (O), 
as  depicted  in  Eqs.  (16)  -  (19)  and  beyond,  are  just  as  beneficial  with  respect  to  Eq.  (23) 
as  they  are  with  respect  to  Eq.  (14). 
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